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ABSTRACT 



Context. The well known radiogalaxy Cen A has been recently detected as a source of very high energy (VHE) y-rays by the HESS 

experiment just before Fermi/LAT detected it at high energies (HE). The detection, together with that of M87, established radiogalaxies 

as VHE y-ray emitters. 

Aims. The aim of this work is to present a lepto-hadronic model for the VHE emission from the relativistic jets in FR 1 radiogalaxies. 

Methods. We consider that protons and electrons are accelerated in a compact region near the base of the jet, and they cool emitting 

multi wavelength radiation as propagating along the jet. The proton and electron distributions are obtained through steady-state 

transport equation taking into account acceleration, radiative and non-radiative cooling processes, as well as particle transport by 

convection. 

Results. Considering the effects of photon absorption at different wavelengths, we calculate the radiation emitted by the primary 

protons and electrons, as well as the contribution of secondaries particles (e*, tts and //s ). The expected high-energy neutrino signal 

is also obtained and the possibility of detections with KM3NeT and IceCube is discussed. 

Conclusions. The spectral energy distribution obtained in our model with an appropriate set of parameters for an extended emission 

zone can account for much of the observed spectra for both AGNs. 

Key words. Galaxies: radio galaxies: individual: Cen A, M87; heavy jet: lepto-hadronic high energy emission 



1. Introduction 

The non-thermal high-energy emission from Active Galactic 
Nuclei (AGNs) has been widely studied in recent years at differ- 
ent wavelength ranges through both satellite-borne and ground- 
based detectors (e.g. Urry & Padovani 1995). Several theoret- 
ical models have been proposed to explain the electromagnetic 
emission of these objects. It is commonly accepted that the high- 
energy radiation is emitted by particles accelerated in the rela- 
tivistic jets launched from the inner parts of an accretion disk 
that surrounds the central black hole. In general, the high-energy 
l" ■ spectral energy distribution (SED) of AGNs presents two char- 
. !^ ' acteristic bumps. The lower energy bump, located at optical to 
S^ X-ray energies, is usually explained as synchrotron emission 
;_( of electrons while the origin of the high-energy peak in the 
Ci ^ SED is still under debate (see e.g. (e.g. Bottcher 2007) for a 
review). Leptonic models attribute this component to inverse- 
Compton up-scattering off synchrotron or external photons from 
the disk and/or radiation reprocessed in nearby clouds (see e.g. 
Katarzyriski et al. 2001; Lenain et al. 2008; Ghisellini et al. 
2005). In hadronic models, interactions of highly relativistic pro- 
tons in the jet with ambient matter and photon fields, proton- 
induced cascades, or synchrotron radiation of protons, are re- 
sponsible for the high-energy photons (see e.g. Miicke et al. 
2000; Mucke & Protheroe 2001; Reimer et al. 2004; Costamante 
et al. 2008; Hardcastle et al. 2009; Romero et al. 1996; Orellana 
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& Romero 2009). There also exist models which are not based 
on the emission of accelerated particles in the relativistic jet and 
assume the production of TeV y-rays in a pulsar-like cascade 
mechanism in the magnetosphere of the black hole (e.g. Neronov 
& Aharonian 2007; Rieger & Aharonian 2008b). 



The recently reported detection by HESS of the nearby ra- 
diogalaxy Cen A (Aharonian et al. 2009) is of great relevance 
since it establishes radiogalaxies as VHE y-ray emitters. Cen A 
is the second non-blazar AGN discovered at VHE, after the 
HEGRA detection of y-rays from M 87 (Aharonian et al. 2003) 
and the later confirmation by HESS (Aharonian et al. 2006). A 
great variety of leptonic and hadronic models has been already 
applied to this kind of sources and a full review is beyond the 
scope of this work. During the first year of operation, Fermi/LAI 
has detected HE emission from Centaurus A (Abdo et al. 2009a) 
and M 87 (Abdo et al. 2009c), providing new constraints to the 
models. 



In this work we present a lepto-hadronic model for the emis- 
sion from FR I radiogalaxies. Section 2 contains a brief descrip- 
tion of observational facts on this type of sources. In Section 3 
we present the outline of our scenario, describing its most rel- 
evant characteristics. Section 4 is devoted to the description of 
the model. In Section 5 we present the application to Centaurus 
A, whereas in Section 6 the results for M87 are given. Finally, in 
Section 7 we discuss the model implications and perspectives. 
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2. FR I radiogalaxies 

According to the unification model of AGNs (e.g. Urry & 
Padovani 1995) FR I radiogalaxies, with their jet axis at a large 
angle with the line-of-sight, are the parent population of BL Lac 
objects whose jets are closely aligned to the line of sight. We 
concentrate here on the only two of them observed until now in 
the VHE range. 

2.1. CenA 

CenA is the closest FRI radiogalaxy (< 4Mpc, Harris et al. 
1984; Hui et al. 1993) and its proximity makes it uniquely ob- 
servable among such objects, eventhough its bolometric lumi- 
nosity is not high as compared to other AGNs. It is very active 
at radio wavelengths presenting a rich jet structure. We can dis- 
tinguish in its structure two components: inner jets at a kpc scale 
and giants lobes covering 10° in the sky. A detailed description 
of the radio morphology can be found in Meier et al. (1989). 
The inner kpc jet has also been detected in X-rays (Kraft et al. 
2002) with an structure of knots and diff'use emission. Recently 
Croston et al. (2009) reported the detection of non-thermal X- 
ray emission from the shock of the southwest inner radio lobe 
from deep Chandra observations. 

The supermassive black hole at the center of the active 
galaxy has an estimated mass of about lO^M© (Neumayer et al. 
2007; Cappellari et al. 2009; Israel 1998) to IO^Mq (Silge et al. 
2005; Marconi et al. 2001). The black hole host galaxy is an el- 
liptical one (NGC5128) with a twisted disk which obscures the 
central engine at optical wavelengths. 

CenA was observed by the Compton Gamma Ray 
Observatory (CGRO) with all its instruments from MeV to GeV 
energies (Gehrels & Cheung 1992; Kinzer et al. 1995; Paciesas 
et al. 1993; Steinle et al. 1998; Thompson et al. 1995). In this 
period, this source exhibited X-ray variability (Bond et al. 1996) 
and also some soft y-ray variability (Bond et al. 1996; Kinzer 
et al. 1995; Steinle et al. 1998). However, Sreekumar et al. 
(1999) found that the EGRET flux was stable during the whole 
period of CGRO observations. 

In 1999 the new Chandra X-ray Observatory took images of 
Cen A with an unprecedented resolution. More than 200 X-ray 
point sources were identified in those images (Kraft et al. 2001). 

Cen A as a possible source of UHE cosmic rays was early 
proposed by (Romero et al. 1996). Recently, the Pierre Auger 
Collaboration reported the existence of anisotropy on the arrival 
directions of UHE cosmic rays (Abraham et al. 2007), remarking 
that at least 2 of this events can be correlated with the Cen A 
position (3° circle). Further works have claimed that there are 
several events that can be associated with Cen A and its big radio 
lobes (Gorbunov et al. 2008; Fargion 2008; Wibig & Wolfendale 
2007) but this correlation is still statistically weak. 

Finally, Fermi/LAT has detected CenA in the first three 
months of survey with a significance above lOcr (Abdo et al. 
2009b). 

2.2. M87 

The giant radiogalaxy M87 is located at 16.7 Mpc within the 
Virgo cluster (Macri et al. 1999). It presents a one-sided jet 
which is inclined with respect to the line of sight an angle be- 
tween 20° - 40° (Biretta et al. 1995, 1999). In addition to its 
bright and well resolved jet, M87 harbors a very massive black 
hole (6.0 ± 0.5) x 10^ Mq (Gebhardt & Thomas 2009) which 
is thought to power the relativistic outflow. Given its proxim- 



ity, the substructures inside the jet could be resolved in the X- 
ray, optical, and radio wavebands (Wilson & Yang 2002). High 
frequency VLBI observations have resolved the inner jet up to 
about 70 Schwarzschild radii (Junor et al. 1999). Along the jet, 
nearly stationary components (Marscher et al. 2008) and fea- 
tures moving at superluminal speeds (Ly et al. 2007; Kovalev 
et al. 2007) were observed (100 pc-scale). 

M87 is also a well-known VHE y-rays emitter (Aharonian 
et al. 2003, 2006; Acciari et al. 2008; Albert et al. 2008) show- 
ing a y-ray flux variability on short time scales with flaring phe- 
nomena in VHE, radio, and X-ray wavebands simultaneously 
(Acciari et al. 2009). Recently, it was detected by Ferm//LAT 
with a significance greater than lOcr in 10 months of observa- 
tions (Abdo et al. 2009c). 

Rapid variability constrains the emission region extent to less 
than a! 56R5, where 6 is the relativistic Doppler factor. Some 
suggested explanations for the VHE y-ray emission were ruled 
out (e.g. dark matter annihilation (Baltz et al. 2000)) At the 
same time, various VHE y-ray jet emission models were pro- 
posed: leptonic (Georganopoulos et al. 2005; Lenain et al. 2008) 
and hadronic (Reimer et al. 2004) ones. However, the location 
of the emission region is still unknown. The nucleus (Neronov 
& Aharonian 2007; Rieger & Aharonian 2008a) , the inner jet 
(Tavecchio & Ghisellini 2008), or larger structures in the jet such 
as the knot HST-1, have been discussed as possible sites of par- 
ticle acceleration (Cheung et al. 2007). 

3. Basic Scenario 

We assume that a population of relativistic particles can be ac- 
celerated to very high energies close to the base of the AGN jet. 
These primary electrons and protons carry a fraction of the to- 
tal kinetic power of the jet L '"', and as they are dragged along 
with the jet, they cool giving rise to electromagnetic emission 
and neutrinos. 

Assuming that a fraction q^ of the Eddington luminosity is 
carried by the jet and a counter jet, the jet kinetic power is 



^(kin) ^ 



qj AnGMhhmpC 



cr-Y 



(1) 



This power can be very high if the jet is launched by a dissipa- 
tionless accretion disk (Bogovalov & Kelner 2010). 

Most of the jet content is in the form of a thermal plasma 
with a constant bulk Lorentz factor Tb. This plasma is initially in 
equipartition with a tangled magnetic field at the Alfven surface 
(zo = 50Rg from the central black hole) (e.g. McKinney 2006). 
The highly disorganized magnetic field has a root mean square 
value B{z) at a distance z from the black hole in the observer 
frame, such that B^{z) := {BI}+{B^,}+{B^}. The magnetic energy 

density for z - Zq is then pmizo) - B^li^Ti) with Bo = B(zo)- 
Equating the magnetic to the kinetic energy density, yields: 



8L' 



(kin) 



\ Hzo)] 



(2) 



Vb 



where za is the distance to the black hole, Vb is the jet velocity, 
and rj(zo) = Zo t^n^j is the radius of the jet assuming that it has 
a conical shape with half-opening angle f|. A widely accepted 
view is that jets are accelerated through the conversion of mag- 
netic energy into kinetic energy (e.g. Komissarov et al. 2007). 
We adopt a phenomenological dependence on the distance to the 



M. M. Reynoso, M. C. Medina, & G. E. Romero: Lepto-hadronic model for FR I radiogalaxies 



black hole for the magnetic field (e.g. Krolik 1999), 



Si,-B,M 



(3) 



Since the density of cold material within the jet decays as z , us- 
ing an exponent me (1, 2) in the above expression implies that, 
as z increases, the magnetic energy decreases more rapidly than 
the kinetic one. The corresponding increase in the bulk Lorentz 
factor is taken into account as described in Appendix A. In the 
following, we will write simply Fb, but it actually depends on z- 
The particle acceleration takes place in a compact but inho- 
mogeneous region of size Az < z-^cc near the base of the jet, at a 
distance z-^cc away from the black hole. The value of z-^cc is fixed 
by requiring the magnetic energy density to be in sub-partition 
with the jet kinetic energy density. This condition enables strong 
shocks to develop (Gaisser 1990). Assuming that at Zacc the mag- 
netic energy is a fraction q^ of its value at zo, it follows that 



zoqn 



(4) 



For example, if Mhh - 10 M©, m - 1.5, and zo - 50Rg, using 
^m = 0.38 yields a distance z-^cc = l32Rg. 

The power injected in the form of relativistic particles (Lrei) 
is considered to be a small fraction (j'rei of the total jet kinetic 
power and the relation between proton and electron powers is 
given by the parameter a such that Lp - a L^ (see e. g. Romero 
& Vila 2008; Vila & Romero 2010). 

These parameters are constrained by the available observa- 
tional data on each source and reasonable theoretical considera- 
tions. In the following section, we describe the procedure used 
to obtain the particle distributions along the jet and the radiative 
output that then arises. 

The existence of heavy jets magnetically driven in AGN 
and microquasars has been supported by several scientific works 
in the last few years. For example, S^dowski & Sikora (2010) 
conclude that mildly relativistic proton-electron jets might be 
formed by magnetocentrifugal launching by inner portion of 
magnetized disks around rotating black holes. This leads to a 
triple-component jet structure: proton-electron component sour- 
rounded by the relativistic pair-dominated sheat. The final speed 
of this centrifugal outflow depends strongly on the disk vertical 
structure (see e.g. Wardle & Koenigl 1993) which is certainly 
unknown. However, for barionic outflows to reach mildly rel- 
ativistic speeds, some inital boost would be neccesary, which 
would be produced by heating or mechanically, by flaring activ- 
ity and/or by radiation pressure. Previous work by Heinz (2008) 
has proved the existence of a relativistic and baryon-loaded jet in 
Cygnus X-1 and that the bulk of the kinetic energy is carried by 
cold protons, as in the case of SS 433. For further information 
about magnetically launched jets see (Spruit 2010), where the 
production, acceleration, collimation, and composition of jets 
are well explained. 



The procedure begins with the calculation of the distribu- 
tion of primary electrons along the jet taking into account syn- 
chrotron and adiabatic cooling. After that, the synchrotron ra- 
diation emitted by the primary electrons can be calculated. To 
check that the electron distribution is consistent with the energy 
loss mechanisms operating, the synchrotron cooling rate must 
be much greater than the inverse Compton (IC) one due to elec- 
trons interacting with the synchrotron photons (SSC). If this is 
the case, it means that the main cooling is due to synchrotron 
radiation, and neglecting the IC energy loss is a valid approx- 
imation to obtain the electron distribution. If the SSC cooling 
can not be neglected, then the transport equation becomes more 
complicated and a different approach is needed (e.g. Schlickeiser 
2009). 

Having obtained the electron distribution, the next step is to 
calculate the distribution of primary protons taking into account 
the cooling due to synchrotron emission, adiabatic expansion, 
pp and py interactions. The latter two types of interactions yield 
the production of secondary pions, muons, and electron-positron 
pairs. These three populations of particles are also described 
with the transport equation, and the radiative output that they 
produce is also considered. According to, e.g. Khangulyan et al. 
(2008) and Pellizza et al. (2010), it can be seen that in the present 
scenario, IC cascades are suppressed by the synchrotron cooling 
of secondary e"^, since the magnetic field is greater than 10 G in 
the regions of the jet where emission takes place. Therefore, we 
neglect the effect of IC cascading and calculate the synchrotron 
emission of the secondary electrons and positrons. 

In this section we present all the relevant expressions used in 
this model. We discuss on the injection of primary particles, the 
relevant cooling rates, the transport equation used, the injection 
of secondary particles, and the emission of photons and neutri- 
nos. 



4. 1 . Injection of primary particles 

At a distance z to the black hole taken in the reference frame 
of the observer, we adopt an injection function in terms of the 
particle energy E' in the co-moving reference frame of the jet: 



Q'i(E\z) :-- 



dNi 



dE'dV'dQ.'dt' 

\2 



(GeV^'s^'cm^^sr"') 



= K, 



/Zaccj ^,-. 



exp 



E' 



p., (max) 



(5) 
(6) 



Here, Ni represents the total number of relativistic electrons 
(/ = e) or protons (/ = p), and the cut-off energy is obtained 
from the balance of particle gains and losses. These processes 
are described below. 

The injection function can be transformed to the observer 
frame by taking into account that 



4. Description of the model 

The model developed for this work is based on the energy distri- 
bution of the different particle populations along the jet. These 
are obtained as solutions of a 1 -dimensional steady-state trans- 
port equation that includes the relevant cooling terms and a con- 
vective one. The radiative output is obtained in the jet reference 
frame, where the particle distributions are isotropic, and the re- 
sult is transformed back to the observer frame. 



E' 



dNi 



p'^ dV'dp'dQ'dt' 



(7) 



is a Lorentz invariant (e.g. Dermer & Schlickeiser 2002; Torres 
& Reimer 201 1). From this, it follows that the particle injection 
in the observer frame is given by 



QiiE,z)^ 



E - mfc' 



2^4 



'2 - m}c^ 



1/2 



Q'i(E',z'), 



(8) 
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where the energy in the jet frame is given in terms of the energy 
in the observer frame E, and the angle 9 between the particle 
momentum p and the bulk velocity of the jet Vb — c/SbZ: 



E' ^tJe- Pb cos e yJE^ - mjcA . 



(9) 



The normalization constant K, is found in each case using the 
power in relativistic species Li^/jI, integrating in the energy, solid 
angle, and volume of the acceleration zone (z e [Zacc, Zacc + Az]): 



Li =1 dE \ dD. \ dV E Qi(E, 

JaE Jatt JhV 



z). 



(10) 



In the jet frame, the cold matter density is 

ndz) 



<{z) = 



Tb 



Relativistic protons in the jet undergo pp collisions with 
these cold protons at a rate 



t'-;„{E',z)^n'^{z)ca-ff\E')Kpp. 






pp 



(16) 



Here the inelasticity coefficient is Kpp ^ 1/2, the corre- 
sponding cross section for inelastic pp interactions can be ap- 
proximated by (Kelner et al. 2006)(Kelner et al. 2009) 



4.2. Accelerating and cooling processes 

To determine the maximum energies, it is necessary to ac- 
count for the particles acceleration and cooling rates t'^ - 
E'-^\dE'/dt'\. 

We assume that in the acceleration zone, particles are accel- 
erated by a diffusive shock-modulated mechanism with a rate 
(e.g. Begelman et al. 1990): 



cr(™i>(£^) = (34.3 + 1.88L + 0.25L2) 



X 10"2^cm^ (17) 



where L = ln(£„/1000 GeV) and £th = 1-2 GeV. 



.-i.p, , ceB'(z) 



(11) 



where < 77 <k 1 is the efficiency of the mechanism. The mag- 
netic field in the jet frame is assumed to be random and with no 
preferred direction. B' represents the root mean square field in 
the jet frame: B'^ = {B'/} + {B'^.}) + {B'/}. As considered also 
in (Heinz & Begelman 2000), the average of each component 
of the magnetic field is supposed to vanish in the jet frame: 
{B'.) = 0. Assuming that there is no electric field in the jet frame, 
the components of the magnetic field in the observer frame are 
related to the ones in the jet frame as B^ = ^hB[,, By = rb^',, 
and B, = B',, which yields' 



B'iz) 



'2r^ + i 



B(z). 



(12) 



where B(z) is given by Eq.(3). 

As for the energy loss processes, particles emit synchrotron 
radiation at a rate 



^ sync(^ '^) 



\m } 



o-jB'^ E' 
nipC Stt mc^ 



(13) 



The lateral expansion of the jet implies an adiabatic cooling rate 
(Bosch-Ramon et al. 2006) in the observer frame given by 



i Z 



(14) 



The density of cold matter at a distance z from the black hole 
is, in observer frame. 



ndz) 



(1 -^rel)»Jj 



m^TTz^ tan^^jVb' 
where the mass loss rate in the jet is 

L*'"'(rb - 1) 



(15) 



' Note that Eq.(12) is consistent with Eq. (3) of (Heinz & Begelman 
2000) if the radial component of their flow velocity can be neglected. 



4.2.1. Synchrotron radiation and Inverse Compton 
interactions 

To consider IC interactions of primary electrons with syn- 
chrotron photons in the jet (SSC), it is necessary to know 
the particle distribution of the synchrotron emitting particles, 
N'(E',z), which is obtained by solving the transport equation 
described below. Actually, the synchrotron emission of electrons 
dominates the low energy photon background which is a good 
target for IC and py interactions. Since in the cases considered 
here, IC cooling is much less efficient than synchrotron cooling, 
the electron distribution can be obtained to a good approximation 
without considering IC cooling. In the jet frame, the background 
radiation density, in units [GeV 'cm"^], has been approximated 
locally as 



«ph(^ph'Z) 



«syn(-E;h'^)rj(z) 



E'u 

ph 



(18) 



where s' is the power per unit energy per unit volume of the 
synchrotron photons. 



^.svn(^Dh'-^) ~ 



^synV^ph 



1 _g-^ssA(£;h,z)' 
TSSA(£;b'^) . 






dE'4nP,ynK(E',z). (19) 



Here, the synchrotron power per unit energy emitted by the elec- 
trons is given by (Blumenthal & Gould 1970): 



P.yn(E;^,E',z) 



^e^B'(z) -gph 
nieC^h E„ 



•JE\ 



lEc, 



d(Ks/3(0, 



(20) 



where K^/^iO is the modified Bessel function of order 5/3 and 



^ _ y/6heB'(z) ( E' 



The factor in parenthesis in Eq. (19) accounts for the effect of 
synchrotron self-absorption (SSA) within the jet, with an optical 
depth 

TSSA(£ph>z)= I ^^'aSSA(£ph-2'(^')), 
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where ff is such that tan 0' - Fb tan 6, and the SSA coefficient is 
given by (Rybicki & Lightman 1979): 



assA 



/ 3 2 r»c 

he r 

2>nE' Jme, 



dE'E"P,UE'E') 



synV^ph> 



dE' 



K(E',z'(l')) 



7/2 



(21) 



After all these considerations, the cooling rate due to IC scat- 
tering for electrons of energy £" in the jet frame can be obtained 
by integrating in the target photon energy E'^ and in energy of 
the scattered photon E'^ (Blumenthal & Gould 1970): 



o 2 4 

t',l(E',z)^- ' 






Vt f 

pi, 

X ' dE;F(q)[E;-E'^^], (22) 



(e.g. Jokipii & Parker 1970; Kirk et al. 1988). The source term 
is given by a function Q'(E',z), which in the case of primary 
particles is given by Eq. (6), while for secondaries it is obtained 
using the parent particle distribution (N'(E',z)) along with the 
secondary particle production rates (see Section 4.3. 1). 

We can solve the transport equation using the method of 
characteristics, i.e., writing 



dz 



dE' 



dN'{E',z) 



b(E',z) Q^^E',z)-[ri^,^'-^]NiE,zy 



(26) 



where the first two terms allow us to find a characteristic curve 
Zc{Ec) for each pair iE',z) of interest. Equating the second and 
third members it follows that 



N'(E',z) 



CdE.^:^ 

Je me,, 

X exp I dE ' 

JE' 
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where £■'*""' is the lowest energy of the available background of 
synchrotron photons, and 

1 (oD^ 

F(q) = 2^1n<7 + (1 + 2^)(1 -q) + -(1 - qyf:^, (23) 



with r; = AE'^^E'limlc^) and 



Ey 



r;£:,(i -£;/£;,) 



-phV 



4.3.1. Pion and muon production 

As high energy protons interact with background matter and ra- 
diation, they produce pions. 

The pion injection due to pp interactions in the jet frame is 
calculated as 



dx _^ _^f t d j^ 



Q'„^p„(K^z)^n'Sz)c r —N'i 
Jo X \ 



E'A r yjEi 
PP 



xFAx,:^\a-ffn^\ (28) 



4.2.2. Proton-photon interactions 

Proton interactions with the photon background can be an im- 
portant cooling process, with a rate appoximated as 



C(E\z)^^J^^dE; 



«;h(^;h'^) 



^ ph 



-2E' 



ph"^/' 



j^^' «'er<'(er)^;r;(er) fr. (24) 



Here, ej^' = 150 MeV and we use the expressions for the cross 

section cr^y and the inelasticity xfy given in Begelman et al. 
(1990). 



4.3. Particle distributions in tine jet 

The energy distribution for each particle population is obtained 
as a solution to the following 1 -D stationary transport equation 
(e.g. Ginzburg & Syrovatskii 1964; Khangulyan et al. 2008): 

dN'(E',z) d[b(E',z)N'(E',z)] N'(E',z) ^,_, , ,. _, 

^b ^ + TTT. + -T^TTT^ ^Q(E ,Z).(25) 



dz 



dE' 



Ta(E') 



This equation includes the effects of convection with a speed 
Vb ~ c, and particle cooling with an energy loss b(E',z) - 
dE' Idt. Particle decay with a timescale Td(E') is also consid- 
ered in the cases of secondary pions and muons. We note that 
the energy E' in the above equation corresponds to the co- 
moving frame, where the particle distributions are isotropic, 
while the spatial coordinate z corresponds to the observer frame 



where 



£: 



' X / \1 +r'x'Hl -X") 



1 -x" 



l-x" l+r'x'^d-x") 



2\l/2 



x|-^+. ^'^^-'""^ .)fl-^l (29) 



£; 



is the distribution of pions produced per pp collision, with x - 
E'JE', B„ = fl'-H0.25, a' = 3.67+0.83L-h0.075L2, / = 2.6/ V^, 
and a = 0.98/ ^/^ (see Kelner et al. 2006, 2009). 

In the same way, the source function for charged pions pro- 
duced by py interactions is 

Q'„„(E',,z)^ f dE'N'^(E',z)aj'fl(E',z) 
Je'„ 

X NAE') 6(E'„ - 0.2E') 
= 5 n;,(5E',, z) w'j;'(5£;, z) NASE'J. (30) 

Here w'^'L' is the py collision frequency defined as (Atoyan & 
Dermer 2003): 



/-*oo 

j'f^(E',z)^:fj L«'£; 



Kh^Kh'^'^ 



^Ph ^,2 
ph 



X2ey„ 
«feV«(e')e', (31) 



py^ 



and the mean number of positive or negative pions is 



(32) 
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This number depends on the probabihties of single pion and unit, per volume unit, per photon energy, and per solid angle. In 
multi-pion production pi and p2 - ^ - Pi- Using the mean in- the observer frame, the corresponding emissivity is given by 



elasticity function Kpy - tj^lupy, the probability pi is 



Qy{Ey,z)^DQ'{E'z\ 



(39) 



P\ 



K2 - Kpj 
K2 - Ki 



(33) where the Doppler factor is D = F^^ (I - /3b cos0) ' and Ey = 



where Ki = 0.2 and K2 = 0.6. 

The injection functions of charged pions given by Eqs. 
(28,30) are used to work out the distribution N'„{E' ,z) by solving 
the transport equation (25). 

For muon injection, we proceed as (Lipari et al. 2007) and we 
consider the production of left handed and right handed muons 
separately, which have different decay spectra: 



DE'y. 

In the case of synchrotron radiation, the emissivity in the jet 
frame is 



yy,syn(-^7'^) 



^nE' 



(40) 



and the IC emissivity is 



dn„-_ 



>t'L 



dE^ 
dn„-. 



X^fij ^n) 



'^R 



dEn 



\*^Li-i ^n) 



rA^ - x) 

E„x(l - r„)^ 

(x - r„) 



®{x-r„) 



E„x(l - r„y 



@(x - r„). 



(34) 
(35) 



Q'y.iciE'z) 



rzc ' " 






dE. 



«;h(£;h'^) 



ph" 



■^ph 



•JE' 



■■Emax N'(E' 7~) 

X I dE' '' ' ' F(q), (41) 

EL.. le 



with X - E^IE„ and r„ - {m^lm„). 

^ The injection function of negative left handed and positive ^^ere r, is the classical electron radius, and we integrate in the 

target photon energy E' ^ and in the electron energies E' between 



right handed muons is 

, dE'J-,!,{E'„ 



dn„ 



^ 5: + !!!££! 

min 2 2 



"f'L, 



v^, 



ph 



x[K-(K,z)^^(E'^;E'^) 

UTlj^+—)ii- 

+N'AK,z)^^(E;;E'„) 



and 

(36) E' = 



F' 



F' 

1 _ y 

£' 



(42) 



(43) 



Given that CP invariance implies that dn„-^^- IdE' — 
dn„.^i_,^JdE'^, and since the total distribution obtained for all With respect to the hadronic contribution to the total emis- 

charged pions is N'„(E'„, z) = A^;, (E'„, z) + N'^ {E'„, z), it follows sion, the photon emissivity due to pp interactions is obtained as 
that the injection for left handed muons is 



, dE'j;},{E'„) 



Q'y^pp(E'y,z)-K(z)cj^^ ^KI^-^,z 



X NUE',,z)^^^(E'^;E'J. (37) 



XFylx, 



E' 



1 ^(inel) 



(T 



PP 



— I, (44) 



And in a similar way, the right handed muons injection is 

Xoo 
, dE:T-},(E'J 

X NUE'„,z)^^^(E'^;E'J. (38) 

Again, the injection functions of muons given by Eqs. 
(37,38) are used to work out the distribution N'(E',z) by solving 
the transport equation (25). 

4.4. Electromagnetic radiation 

The main radiative processes considered in this work are: syn- 
chrotron radiation, IC emission, pp and py interactions. For each 
emission process we can calculate the injection of photons or ra- 
diation emissivity, which is given by the general expression of 



where the function Fy{x, E') is the same as defined by Kelner 
et al. (2006), for a proton energy E' = E' /x. 

For the y-ray emissivity due to py interactions we use the 
expression 

Qr^Pr(E'r^^-i^l,§^K(E',z) 

"" ri/^>;h(£;h'^)'I'('7.r'^)> (45) 

ph 

4E' E' 

where E'y = xE', T]py = -3-4^, 



7 (min) 
-ph 



2 4 

m c X 
AE', 



(46) 



and with the function (^\ripy,x\ as tabulated in Kelner & 
Eq. (5) and represents the number of photons produced per time Aharonian (2008). 



M. M. Reynoso, M. C. Medina, & G. E. Romero: Lepto-hadronic model for FR I radiogalaxies 



4.5. Internal and external photon absorption 

We consider the absorption of the radiation due to photoioniza- 
tion processes at eV energies, and also due to e^e^ creation by 
77 and -yN interactions at high energies in the source. 

For yN interactions, we take as target the material along the 
line of sight corresponding to each particular object, which is 
measured through the column density of neutral hydrogen Nh- 
The absorption cross section o-y^/ in this case is taken as in Ryter 
(1996) for Ey < I keV, assuming that the medium is composed 
by atomic hydrogen and dust within galactic abundances (see 
Fig. 1). This cross section includes the effects of photoionizaton 
for Ey > 13.6 eV and scattering with dust below this energy. In 
the present work we do not consider recombination, the inverse 
process of photoionization (see Appendix B ). The resulting op- 
tical depth for jN interactions is approximated as 



TyH(Ey) = NnO-yNiEy). 



(47) 



It can be seen from Fig. 1 that besides the large absorp- 
tion edge corresponding to the ionization energy of hydrogen, 
three aditional absorption edges are included in the cross sec- 
tion. They correspond to different ionization energies of helium, 
K-shell electrons of oxygen, and iron. In the present context, 
as can be seen below, we find no significant features associated 
with these aditional absorption edges (e.g. Cruddace et al. 1974; 
Ghisellini et al. 1999). The possible re-emission of lines is be- 
yond the scope of this work, and could be studied, e.g., as in 
Ghisellini et al. (1999). For energies above 1 keV we take the jN 
cross section from Amsler et al. (2008), to account for Compton 
scattering and e^ production. In Fig. 1 we show the cross section 
used for 7 A^ interactions 
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Fig. 1. Cross section for jN interactions. The absorption edges 
corresponding to the different elements are indicated. 

In the case of 77 interactions, the main radiation target is 
considered to be the synchrotron photons inside the jet, which 
are characterized by a radiation density n'^(E'^^,z)- If a dissi- 
pationless accretion disk is present, no other significant photon 
field is relevant for gamma-ray absorption. The cross section is 
given by 



cryy(E'y,E;^)^^(l-/3l)x 



2;S,(^^2) + (3-/^)ln(i^| 



with 



/3c 



2mic' 



2^4 



£';£'pl,(l - COS Oyy) 



1/2 



(48) 



where Oyy is the angle of interaction between the incident pho- 
tons. For 7-rays produced at a position z in the jet, we assume 
that they undergo collisions over a length r^{z)l s,mff inside the 
jet with a target radiation field that is isotropic in the jet frame. 
The corresponding optical depth is calculated using the general 
expression of Gould & Schreder (1967) to obtain 



rint) 



(£;,z) = 27r y"" dl' \ dx'(l-x') 

^ J 2.y ^Vpi^^K^' z'(l'))cryy(E'y, £;,). (49) 



In Fig. 2 we show the optical depths obtained for the case 
of Cen A and M87, Tjn and Tyy, the latter evaluated within 
the injection zone and also outside it. The absorbing photons 
outside the injection zone are those of synchrotron emission of 
protons and secondary e=^, which in that case are also included 
in n'^(E'^^,z)- It can be seen from this plot that very important 
absorption occurs at low energies by photoionization, and also 
at VHE through internal 77 interactions. The latter effect takes 
place mainly at the injection zone, since for z > Zacc + Az, the 
optical depth is much lower. 

Taking into account the total optical depth Ty{Ey,z) = 
Tyy(Ey,z) + Ty^iEy), thc differential photon flux at Earth, i.e., 
the number of photons with energy Ey, per unit energy, per unit 
area, per unit time, can be calculated as 



d%{Ey) 
dEy 



^ \ dVQy(Ey,Z)tX^[-Ty{Ey,Z)\ (50) 



where d is the distance from the source to Earth. 



4.6. Neutrino emission 

Neutrinos arise from direct pion decays plus muon decays with 
a total emissivity 

Q',/E',z)^Q'„^,^(E',z) + Q'^^,^(E',z), 

which correspond to the observer frame and transforms accord- 
ing to Eq. (39). The contribution from pion decays (n^ — > ju^v^, 
n^ — > yu^v^) is represented in the co-moving frame by 



Q'(E',z)= dE„T-„\(E„)K{,E„,z) 
Je' 



^m^^Jj^^ (51) 



E„{l-r,) 



with X = E' /E„ and T^j d - Jn 2-6 x 10 s 



The neutrino emissivity for muon decays (yu 
— > e^VgV^) can be calculated as 



^e^ji. 



Q',^^^{E',z) = X r ^T-^!/E,)N;_(E„z) 
,_, Je' ^p 



5.24, 
3 3 



(52) 
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Fig. 2. Optical depth for Cen A and M87 at different distances 
form the black hole. 

In this expression, x - E' jE^, ^\2 - f^^'^, Tfi^d = yf, 2.2 x 10"''s, 
and /i3 4 = /ig'^ (see Lipari et al. 2007). The differential neutrino 
flux can be computed using a expression analogous to Eq. (50). 
The synchrotron cooling of fis and tts, which is taken into 
account in our model, will be traduced on a lower neutrino emis- 
sivity as it was previously shown in Reynoso & Romero (2009). 
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Here, the photopair cross section o-^y is considered as in 
(Begelman et al. 1990). 

As discussed above, the inelastic pp and py collisions pro- 
duce pions that decay giving muons. The decay products of 
muons include electrons and positrons. The decay spectrum of 
the electrons and positrons are taken to be equal to that of muon 
neutrinos and anti-neutrinos (see e.g. Gaisser 1990). Thus, the 
injection of electrons and positrons from muon decay after pp 
and py interactions is considered following Schlickeiser (2002). 

As for the e^e^ pairs produced in yy interactions, we con- 
sider the injection according to Aharonian et al. (1983), 
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(55) 



where y^ = E' KnieC^), o) = E'^JinieC^) , and e = EylinieC^). 
The soft photon density n' ^ refers to the synchrotron radiation 
of electrons, while the density of higher energy photons from 
IC, pp, and py interactions, is given by 



n'y(E'y, z) 



An 



Q'yjc(E'rZ) + 6;.„,(£;,z) + Q'y,py(E'y,z) 
c/(ztan^j) + f-i(£' z) 



Here, the annihilation rate of gamma-rays due to pair production 
is given by 



t;l(E'y,z)^2ncf dx'(l-x') 



r*oo 

■J 2^ 



dE'h<h(Kh^z)o-yy(E'E'). (56) 



ph^ 



The distribution of secondary electrons and positrons in the 
jet is found solving the transport equation. This allows us to cal- 
culate the synchrotron and IC emission produced by these parti- 
cles. 



4.7. Secondary e^e production 

We consider the production of secondary electrons and positrons 
in the jet through py, pp, and yy interactions. The inelastic 
py collisions with can produce e^e^ pairs apart from pions as 
discussed above. We consider this direct production of pairs 
through the injection function in the jet frame 



QU^„y(E',z)^2 dE;N;(E;,z) 

Jm,,c^ 



XOJ 



'rl(E',,z)6[E' 



E'\, (53) 



where the ojpy^e*(E') is the collision frequency given by 
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«;h(^;h'^) 
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ph 
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de'cT^plie) e. 



(54) 



5. Application to Centaurus A 

In this section we present the results for Centaurus A, which 
were obtained applying our model with the set of parameters 
listed on Table 1 . The cooling rates for high energy electrons and 
protons for this configuration are shown in Fig. 3, while the ob- 
tained electron and proton distributions N'e(E',z) and N'(E',z) 
are presented in Fig. 4. 

We can see in Fig. 3 that the relevant energy looses are due 
to synchrotron cooling for electrons and py interactions for pro- 
tons. Since the plot corresponds to the injection zone (z = Zacc), 
py interactions are favored by a large density of target photons 
corresponding to the synchrotron emission of electrons. These 
photohadronic interactions are not so important for protons out- 
side the injection zone. 

It is important to notice that the maximum energy achiev- 
able for the protons in this context is of 2 x 10^ GeV. These 
protons can not account for the UHE cosmic rays detected by 
the Pierre Auger Observatory in the direction of Cen A but it 
gives a hint to other explanations. If neutrons of the same en- 
ergy are produced by pion photoproduction inside the jet, these 
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Fig. 3. Accelerating and cooling rates for electrons {left) and for protons (right) at a distance z = Zacc from the central engine of Cen 
A. In the case of electrons, the light-grey dashed line correspond to the first calculation of the IC cooling rate taken into account only 
synchrotron cooling to obtain the electron distribution. A second approximation to the electron distribution leads to the IC cooling 
rate indicated in white circles. 



Table 1. Model parameters for Cen A 



Parameter 



Value 



Mbh: black hole mass 

Rg'. gravitational radius 

L . jet kinetic power at zo 

g'ji ratio 2L*'"ViEdd 

FbCzo): bulk Lorentz factor of the jet at zo 

d: viewing angle 

^j! jet's half-opening angle 

(/reP jet's content of relativistic particles 

a: hadron-to-lepton power ratio 

Zo'- jet's launching point 

q^'- magnetic to kinetic energy ratio at Zacc 

Zacc^ injection point 

Az: size of injection zone 

S(Zacc): magnetic field at Zacc 

m: index for magnetic field dependence on z 

s: injection spectral index 

7/: acceleration efficiency 

Ep : mmimum proton energy 

Ef : minimum electron energy 

£(max) ■ - , 

p : maximum primary proton energy 

Ef' : maximum primary electron energy 

Nh'- column dust density 

ndZacc)'- cold matter density inside the jet at z-m 



10** Mo 

1.47 X 10"cm 

6.28 X 10'*''erg s"' 

0.1 

3 

30° 

2.5° 

0.05 

0.025 

50 Rg 

0.38 

132 Rg 

2 Zacc tan fj = 11.5^g 
3065 G 

1.5 
1.8 
10-2 

3 GeV 
0.1 GeV 

2 X 10^ GeV 
200 GeV 
10^3 cm-2 

3 X 10** cm-3 



could been beamed along the jet and decay into protons near the 
outer radio lobes, where they would be re-accelerated up to the 
observed energies by Auger. Other possibilities were proposed 
by Rieger & Aharonian (2009), with particles that go through 
shear acceleration inside the jet, and by Romero et al. (1996), 
who considered that the UHE cosmic rays are produced at the 
extended radio lobes. If this is the case, it has been estimated by 
Hardcastle et al. (2009) that the minimum jet power required is 
roughly L > 3.5 x lO'^^ergs"', which is consistent with the 
value assumed here. 

Fig. 4 shows that the electron distribution drops quickly with 
z at the end of the injection zone, while the proton distribution is 
also important further along the jet. This is a consequence of the 
different energy loss rates of the particles: electrons cool very 
rapidly emitting synchrotron radiation, and protons lose energy 
at a much lower rate as they propagate along the jet. 

The SED of Cen A (see Fig. 5) includes the HESS spectrum 
in the VHE range together with data from CG/JO/COMPTEL 
(Steinle et al. 1998), RXTE and INTEGRAL (Rothschild et al. 
2006). We also include data from HST/mCMOS and WFPC2 
(Marconi et al. 2000), SCUBA at 800 iJ.m from (Hawarden et al. 
1993), ISO and SCUBA (450 ^um and 850 //m) (Mirabel et al. 
1999), XMM-Newton, Chandra spectra and the Suzaku data pre- 
sented in (Abdo et al. 2010). The Fermi/LAI correspond to 
Cen A core as given in (Abdo et al. 2010), and, finally, we also 
show the old data from the NASA Extragalactic Database. 

The Ferm//LAT analysis revealed that the high-energy spec- 
trum is non-variable over the first ten months of scientific opera- 
tion of the instrument. This steady behavior is also supported by 
the HESS experiment, which also reported a constant flux from 
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Fig. 4. Distributions of primary particles as function of the energy and the distance to the core inside the CenA jet.Le/f: elec- 
trons; right: protons. 




CenA even though CGRO data shows some variability with at 
least two emission states during the period 1991-2000. 

Since the set of measurements composing the SED of Cen A 
is quite inhomogeneous in time and angular resolution, so one 
should be very careful in attempting to interpret any fit of all the 
spectrum simultaneously. In particular, the data that define the 
bump in the hard-X-rays (~0. 1 Me V) have been taken more than 
10 years ago with a poor angular resolution and long integration 
times. The lack of a good spatial resolution makes impossible to 
distinguish the emission components (jet, nucleus or other radi- 
ation sources). Furthermore, we note a discrepancy between the 
flux normalization of FermifLAI and HESS, which is not yet 
fully understood. 

Nevertheless, we have obtained using our model a spectral 
energy distribution which is basically consistent with the multi 
wavelength emission from Cen A. 

An important role is played by absorption due to photoion- 
ization interactions in the surrounding dust. This arises due to 
the large value adopted for the column density of neutral hy- 



drogen, Nh = 10 cm" (see, e.g. Evans et al. 2004; Morganti 
et al. 2008). A drastic modulation is then imprinted in the elec- 
tron synchrotron spectrum, which is responsible for the whole 
emission in the broadband range 10"^ - lO^eV. Such situation 
is possible if electrons can be efficiently accelerated {rj = 0.01) 
to high energies with a rather flat spectral index {s - 1.8). The 
internal -yy absorption does not modify significantly the pp con- 
tribution of gamma-rays, since it is only important within the 
injection zone, which is a negligible region of the jet compared 
to the one in which pp collisions occur (see Fig. 4). 

We evaluate the accompanying neutrino output using the 
same set of parameters. The obtained differential flux, weighted 
by the squared energy is plotted in Fig. 5 together with the es- 
timated sensitivity of KM3Net for one year of operation (Katz 
et al. 2009). We can see in this figure that the neutrino signal 
produced by pp and pj interactions would not be observable 
by KM3NeT detector in just one year of observation. It is to 
be noticed, however, that the sensitivity level shown in the plot 
actually corresponds to a neutrino spectrum with an E^^ depen- 
dence, which is a bit steeper than our case. Hence, we can expect 
that the actual sensitivity for our flux will be better than what is 
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Fig. 6. Differential neutrino flux weighted by the squared energy 
as predicted by the model for Cen A. Dotted line: Approximate 
KM3Net sensitivity for 1 year operation. 



shown here, so with two or three years of data taking it would be 
possible to achieve detection. 



6. Application to I\/I87 

In the case of M87, we have used the set of parameters listed on 
Table 2 as well as the multi wavelength data available. We show 
the cooling rates for high energy electrons and protons in Fig. 7, 
and the electron and proton distributions N'^{E',z) and N'(E',z) 
in Fig. 8. 

We can see in Fig. 7 that the relevant energy looses are due to 
synchrotron cooling of electrons and protons. Protons also cool 
significantly through adiabatic losses, and although pp interac- 
tions are not dominant, they can give an important radiative out- 
put. From the balance between this radiative loss rate and the ac- 
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The rest of the data points correspond to the references given in the text. 



celeration rate, the maximum energy achievable for the protons 
inside the jet is about 9 x lO*" GeV. These protons do not escape 
from the source. The only possibility to obtain a detectable cos- 
mic rays flux at the earth coming from M87 should depend on 
neutrons produced by pion photoproduction processes and rela- 
tivistically beamed along the jet direction and Doppler boosted 
in energy. 

As in the case of Cen A, Fig. 8 shows that electrons radiate 
only at the base of the jet loosing all their energy rapidly while 
protons survive longer inside the jet, emmiting enough radiation 
to account for the VHE y observed emission. 

In Fig 9 we plot the contributions to the SED of M87 ob- 
tained with the different emitting proceeses in our model. We 
also show the corresponding multi wavelength observational 
data. The HESS spectrum included in the VHE range corre- 
sponds to a low state from the 2004 observing season (Aharonian 
et al. 2006). Data from Fermi/LAT (Abdo et al. 2009c), 2009 
MOJAVE VLBA 15 GHz (Kellermann et al. 2004) and Chandra 
X-ray measurements of the core (Harris et al. 2009) are also 
shown. We also include the 3cr upper limits of the integrated 
emission in three hard X-ray bands based on the Swift/BAT 
dataset (Ajello et al. 2009). Historical measurements of the core 
are also plotted: VLA 1.5, 5, 15 GHz (Biretta et al. 1991), 
IRAM 89 GHz (Despringre et al. 1996), SMA 230 GHz (Tan 
et al. 2008), Spitzer 70, 24 m (Shi et al. 2007), Gemini 10.8 m 
(Perlman et al. 2001) and HST optical/UV (Sparks et al. 1996). 



We have obtained, using our model, a spectral energy distri- 
bution which is basically consistent with the multi wavelength 
emission from M87. The spectrum detected by the Fermi/LAT 
collaboration connects smoothly with the low-state TeV spec- 
trum detected by HESS. Our model reproduces this connec- 
tion effectively with the y-ray emission produced by pp colli- 
sions. Contrary to the case of Cen A, this radiative output for 
M87 is possible if particles are not so efficiently accelerated 
{rj - 1.5 X 10"^) to high energies with a rather steep spectral 
index {s - 2.4). 

The flaring state, which is not shown in Fig. 9, can still be 
reproduced in our model by y-rays from pp interactions if we 
assume a flatter spectrum of injected protons and a negligible 
power in relativistic electrons. 

The accompanying high energy neutrino flux was calculated 
using also the set of parameters of Table 2. We show the differ- 
ential neutrino flux weighted by the squared energy as a function 
of energy in Fig. 10, where we also plot the estimated sensitivity 
of IceCube for one year of operation (Ahrens et al. 2004). 

As can be seen been from Fig. 10, the neutrino signal asso- 
ciated to M87 would be hard to detect. However, the sensitivity 
shown corresponds to a neutrino spectrum with an E^"^ depen- 
dence, i.e., not so steep as our predicted one. Hence, a more de- 
tailed study could help to assess the real detectability of this type 
of signal (Neronov & Ribordy 2009b, a). 
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7. Conclusions 

We have presented here a lepto-hadronic model describing the 
particle propagation and interactions in the inclined jets of ra- 
diogalaxies. The scenario includes an acceleration zone placed 
near the jet base, where relativistic particles are injected. The lo- 
cation of this injection zone is fixed in the jet at a distance z-^cc 
from the black hole by requiring that the magnetic energy den- 
sity is less prominent than the bulk kinetic density. This favors 
shock formation. 



One very important aspect of the present model is that it does 
not consist of a one-zone treatment, since particles are allowed 
to convect along the jet, away from the injection zone. To ob- 
tain the particle distributions, we have used a stationary one- 
dimensional transport equation that accounts for energy losses, 
escape, and convection of particles. The resulting distribution 
of primary electrons is important mainly in the injection zone, 
since they undergo a rapid synchrotron cooling. Instead, pro- 
tons are convected away along the jet, loosing energy gradu- 
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ally (mainly by adiabatic cooling, synchrotron radiation, and pp 
interactions). Also they undergo significant py collisions with 



the electron-synchrotron photons at the injection zone. The sec- 
ondary pions and muons produced in the hadronic interactions 
are also considered, and their distributions are found using the 
mentioned type of transport equation. This allows to account for 
the synchrotron cooling suffered before the decay of this sec- 
ondary, transient particles. 

We have found that the radiative output obtained can account 
for much of the observational data for the two radiogalaxies 
M87 and CenA. As seen in Figs. 5 and 9, the model allows 
for the possibility of explaining the VHE emission by means 
of PP interactions. The hard X-ray data corresponds to elec- 
tron synchrotron radiation whereas the soft-y emission is well 
reproduced by proton synchrotron and IC emission. Previous 
hadronic models (e.g. the synchrotron proton blazar of Miicke 
et al. (2003)) could also fit the Fermi energy range for M87 
and the hard spectrum measured by HESS (consistent with 
VERITAS and MAGIC measurements), according to what it is 
mentioned in Abdo et al. (2009c). A previous one-zone lepto- 
hadronic model for Cen A was proposed by Orellana & Romero 
(2009) just before HESS and Fermi/LAT data was available. 

In the present work, we have included the effects of ab- 
sorption due to photoionization interactions in the surrounding 
medium. This causes a drastic modulation in the electron syn- 
chrotron spectrum, especially for CenA, which is responsible 
for the whole emission in the broadband range ~ 10"^ - lO^eV. 
In the case of M87, this type of absorption has a minor impact, 
since the column density is much less (Nh - 2 x 10^'^cm"^, 
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Table 2. Model parameters for M87. 



Parameter 



Value 



Mbh: black hole mass 
Rg. gravitational radius 

L*"': jet power at zo 

qf. ratio 2L*'"VLEdd 

T],(zo): bulk Lorentz factor of the jet at zo 

9: viewing angle 

^j! jet's half-opening angle 

(/jei: jet's content of relativistic particles 

a: hadron-to-lepton power ratio 

qtn'. magnetic to kinetic energy ratio at z^cc 

Zq: jet's launching point 

Zacc: injection point 

Az: size of injection zone 

S(zacc): magnetic field at Zacc 

m: index for magnetic field dependence on z 

s: injection spectral index 

7/: acceleration efficiency 

£(min) • • , 

p . mmimum proton energy 

Eg . minimum electron energy 

£(max) - - , 

p : maximum primary proton energy 

E^"^ : maximum primary electron energy 

Nh'- column dust density 

Wcfeacc): cold matter density inside at Zacc 



6 X 10" Mo 

8.87 X 10''*cm 
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4 

30° 
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40 
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143 i?g 

Zacctanfj/4=11.5i?g 

407 G 

1.5 

2.4 

1.5 X 10^^ 

3GeV 

65MeV 

9 X 10* GeV 

25GeV 

2 X 102"cm-2 

3.9 X 10* cm-3 



Lieu et al. (1996)). The internal absorption of gamma-rays is 
also taken into account, being important only within the injec- 
tion zone, where the synchrotron emission of primary electrons 
provide an important absorbing target. 

It is important to note that the same model can be used to 
describe different types of jet. In terms of the broadband photon 
emission, M87 data shows a similar luminosity level for high 
and the low energy ranges. Since the VHE emission is basi- 
cally determined by pp interactions, we need a large power to 
be injected in relativistic protons as compared to that injected in 
electrons. This is because electrons cool completely, mainly by 
synchrotron emission, and protons undergo an important adia- 
batic cooling. Hence, in the case of M87, we need a high value 
for the proton-to-electron power ratio: a - 40. A different sit- 
uation arises for Cen A, where the highest luminosities corre- 
spond to energies Ey < lO^eV, and the HE and VHE luminosi- 
ties are comparatively lower In this case, we need more power 
to be injected in electrons than in protons (a - 2.5 x 10"^) to 
allow the broadband electron-synchrotron spectrum to approxi- 
mately account for all the data for Ey < lO^eV. Then, for Cen A 
a lower luminosity in protons is enough to reach the level of 
the observed HE and VHE radiation through proton-synchrotron 
and pp interactions, respectively. A further difference between 
Cen A and M87 is in the slope of the injected particle distribu- 
tions: for Cen A, we have a quite flat injection {s - 1.8), while 
for M87 we find a steeper injection {s - 2.4). Also the accelera- 
tion efficiency is greater for Cen A than for M87. This also has an 
impact in the neutrino spectrum produced, which for Ey > ITeV 
appears more difficult to be detected for M87 than for Cen A. 

Finally, we remark that a possible improvement of the 
present treatment, which is left for future work, should be the 
inclusion of a time-dependent injection to account for flaring 
states. 
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Appendix A: Bulk Lorentz factor evolution 

The sum of the kinetic plus magnetic energy in the jet is 



dznz^X^ \Pkiz) +Pm(z)] = 2p,^(zQ)dznzl)(^ 
where ;tf = tan^j. Then, 
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and given that 
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Squaring, we obtain the equation 

rtOW)^^ 

Tb + 1 
which has the analytic solution 



(A.6) 
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(A.7) 



For Cen A and M87 we obtain the results shown in Fig. A.l. 

Appendix B: Absorption by photoionization in the 
surrounding medium 

Here we discuss on the absorption of photons once they espape 
from the jet and propagate in the ambient medium filled with hy- 
drogen gas. In principle, photons with energies above the Lyman 
edge, 13.6 eV, can ionize the neutral hydrogen atoms, and the 
recombination of electrons and protons can also occur It is ex- 
pected that the vicinity of the jet will be highly ionized, but since 
the emitted radiation can not ionize an infinite volume of gas, 
then the gas must become neutral at some point (see e.g. Spitzer 
1978). 

A crude estimation of the size of the ionized zone can be 
attempted assuming a spherical symmetry of the emitting zone 
and of the radiation emitted. Considering that the gas extends 
upto a distance Rh from the emitting zone, we adopt an uniform 
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Fig. A.l. Bulk lorentz factor of the jet as a function of the dis- 
tance to the black hole obtained for Cen A (left panel) and M87 
(right panel). 



density of hydrogen rin - Nh/Rh- This is actually the density of 
ionized hydrogen plus that of neutral hydrogen: 



riH ^np + riHi- 



(B.l) 



The size of the ionization zone can then be estimated through 
the so-called Stromgren radius rs (e.g. Stromgren 1939), which 
is defined by equating the rate of ionizing photons Ny to the num- 
ber of recombinations per unit time: 



-/Vy — fffec 
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47T 



(B.2) 



Here, the number density of free electrons is taken to be Wg ? 
within the ionized zone, and the rate of ionizing photons is 



nn 
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fdV f 

Jv Je-, 



dEAjiQy^^-s 
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which in our cases corresponds mainly to synchrotron radiation 
of electrons. The recombination coeficient can be taken as a\^l^ « 
2.6xl0"-'cm-'s"' for r = lO'^K and neglecting recombinations to 
the energy level n - \, since this state is ionized again on the spot 
(e.g. Krolik 1999; Hunt & Hirashita 2009)). It can be checked 
that in the cases studied here, the Stromgren radius is less than 10 
pc, which is assumed to be much less than Rh ■ If the ionization 
fraction is x = rip/tiH, then the density of neutral hydrogen is 
«//(! - x), which is an increasing function of the distance within 
the ionized region. The optical depth at a distance R from the 
source and along the photon path can be integrated as 



TyN(Ey,R) 



dly nnil - x)o-yN(Ey) 



(B.3) 



According to Petrosian et al. (1972), the fraction of neutral 
gas reaches unity at the Stromgren radius or even less in the 
prensence of dust. If the rs » Rh we can approximate the opti- 
cal depth at Rh as 



TyN(Ey) X 0-yN(Ey)NH, 

which is the expression we use in this work. 
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